Deuteron A^A^*(1440) components from a chiral quark model 
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We present a nonrelativistic coupled-channel calculation of the deuteron structure including AA 
and A'^A^*(1440) channels, besides the standard A^'A'^ S and D-wave components. All the necessary 
building blocks to perform the calculation have been obtained from the same underlying quark 
model. The calculated NN*{14:40) probabilities find support in the explanation given to different 
deuteron reactions. 
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Since the discovery of nucleon structure, nucleon reso- 
nances have attracted considerable attention from theo- 
rists and experimentalists. An important effort was made 
to choose and design experiments that could probe the 
presence of resonance configurations in different nuclear 
systems. The basic idea lies in the observation that a 
small fraction of the nucleons will be internally excited 
and therefore present as virtual resonances in every nu- 
cleus. This may happen even at low energies due to the 
possibility of exciting internal nucleon degrees of freedom 
according to the process NN NN* (or TV* TV*) 
NN involving intermediate iV*'s. As a consequence, the 
many nucleon wave function should be supplemented by 
configurations involving one or several nucleons in an ex- 
cited baryon resonance state. 

If the virtual iV*'s exist in bound nuclear states, one 
expects them to play an important role already in the 
bound two- nucleon system, the deuteron Q. The most 
prominent low-lying even-parity nucleon resonances are 
the P33, the A(1232), and the Pn, the N*{1U0) or Roper 
resonance. The 7V*'s contribute predominantly to the 
nucleon short-range correlations and enhance the high 
momentum components of the nuclear two-particle den- 
sity. Being the deuteron isoscalar, the energetically low- 
est state A^A is forbidden, and therefore the AA and 
]VjV*(1440) components would be the relevant nonnucle- 
onic configurations. The admixture probabilities of these 
exotic states are small due to the low nuclear density 
and the rather high resonance-nucleon mass difference. 
Nonetheless, they have been advocated long ago to im- 
derstand elastic proton-deuteron backward scattering at 
energies above pion threshold or the angular distri- 
bution of deuteron photodisintegration at energies above 
Ej = 100 MeV 1^. Recent calculations have renewed 
the interest on these nonnucleonic components as they 
could be indirectly observed in several reactions as for 
example antiproton-deuteron annihilation Q , subthresh- 
old antiproton production Q or ^ dp processes 
Although the evidence for resonance configurations in the 
deuteron from such processes is indirect, is suggestive and 
encouraging. 

The treatment of nucleon resonances in the nuclear 
wave function can be done in different ways. One pos- 
sibility has consisted on keeping only nucleons (and no 
resonances) in the nuclear wave function and using effec- 
tive operators. However, it would be surprising if such an 



ad hoc procedure would phenomenologically account for 
any detailed internal dynamics, because in the hard-core 
region two nucleons are likely to excite each other and 
thus mutually probe their internal degrees of freedom. If 
one wants to turn the attention to the virtual contribu- 
tion from these nucleon excited states, one has to include 
explicitly N* transition potentials in a coupled channel 
calculation. 

When performing a coupled channel calculation for the 
deuteron based on effective baryon-baryon potentials, a 
problem immediately arises. If one uses for the nucleon- 
nucleon (NN) channel an effective potential which is 
fitted to the NN scattering data, it will already in- 
clude contributions from intermediate N*'s, and thus 
one would obtain too much attraction at medium range. 
Therefore one has to modify the normal nucleon-nucleon 
potential and weaken the intermediate range attraction 
in order to account for the additional attraction from the 
explicit dispersion contribution to the potential with in- 
termediate N*'s. Such a procedure usually introduces an 
unwanted model dependence on the results obtained. 

There are multiple examples in the literature of these 
type of calculations. Haapakoski and Saarela [Q stud- 
ied AA components on the deuteron changing in an ad 
hoc manner the intermediate range attraction of the cen- 
tral Reid soft-core potential until the deuteron binding 
energy was fit to the experimental value. The tensor 
force was not changed. Arenhovel, Danos and Williams 
Q studied A^iV*(I440) configurations using perturbation 
theory with the one-pion-exchange (OPE) potential and 
Hulthen NN deuteron wave functions. The singular na- 
ture of the potentials required a cutoff factor and the 
results were found to be rather sensitive to the cutoff 
chosen. Weber |^ and Nath and Weber calculated 
NN* components by means of an OPE potential in mo- 
mentum space. The resulting probabilities were again 
found to be rather dependent on the high momentum 
suppression factor. Finally, Rost performed a cal- 
culation where the existence of resonance components is 
taken into account by a modification of the intermediate 
attraction of the Reid hard-core potential. As a conse- 
quence, the study of nonnucleonic configurations in the 
deuteron based on standard meson-exchange NN poten- 
tials depends on two basic assumptions: on one hand, 
the hypothesis done to modify the intermediate range 
attraction of the NN interaction, on the other hand, the 
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specific transition potential to the resonance configura- 
tions used. 

During the last decade an important effort has been 
devoted to the qualitative and quantitative understand- 
ing of the nucleon-nucleon interaction based on quark 
degrees of freedom. The first attempts of the early 80's 
based on the one-gluon-exchange (OGE) force Ina 
evolved to the so-called hybrid quark models Fsf - 
The original OGE was supplemented by Goldstone-boson 
exchanges at the level of quarks and, thus, fully quark- 
model based NN potentials were obtained. The OGE 
provided the first explanation of the repulsive core of 
the NN interaction. Although it is known that quark- 
antisymmetry on the Goldstone-boson exchanges induces 
the same effect the OGE still results crucial to reg- 
ularize the short-range part of the Goldstone-boson ex- 
change tensor force. The pseudoscalar-boson exchange 
tensor coupling is fundamental to reproduce the NN 
phase shifts, and, finally, the scalar-boson exchange pro- 
vides the necessary medium-range attraction. The two- 
baryon problem based on quark-quark interactions is 
solved by means of the resonating group method (RGM) 
in a coupled-channel scheme considering usually NN, 
A A and hidden color- hidden color components. The im- 
portance of particular configurations, like A'^A, to de- 
scribe the low orbital angular momentum partial waves 
has been recently demonstrated Such calculations 

present the great advantage that any baryon-baryon in- 
teraction can be determined, once the NN potential 
has been fixed, in a completely parameter-free way. 
Therefore, the quark-model framework provides an ad- 
equate scheme to study nonnucleonic configurations on 
the deuteron without the aforementioned uncertainties 
appearing in meson-exchange models. 

In the present work we want to focus on the influ- 
ence of the most important nonnucleonic channels on the 
deuteron properties. In order to have a consistent calcu- 
lation we will make use of baryonic potentials constructed 
from the same underlying quark-quark interaction po- 
tential so that no parameters are fitted independently 
for the different channels. The quark-model parameters 
have been previously determined to describe the deuteron 
binding energy and S'-wave NN phase shifts with enough 
accuracy so that conclusions on the role played by the dif- 
ferent channels, and thus, by the A* (1440) resonance or 
A can be inferred. 

Within the quark model framework, this problem has 
already been partially undertaken using different approx- 
imations. The possible effects of a bigger Hilbert space 
were considered in Ref. through the formation of 

six-quark bags at short-distances. Ref. proposed an 
alternative formulation in terms of quark-shell configu- 
rations, later on projected onto physical channels. Ref. 

proposed an indirect calculation where the nonnucle- 
onic configurations were not explicitly considered for the 
calculation of the two-body system. As a consequence, 
these methods rely again on several hypothesis that could 
hide physical conclusions. In Ref. pll] the influence of 



A' and A resonances on the A^A^ interaction has been 
studied. The most significant contribution was obtained 
from channels involving A^ and A ground states, although 
a quantitative calculation of the nonnucleonic configura- 
tions was not performed. 

Among the various constituent quark models already 
available in the literature, to have a reliable prediction for 
the nonnucleonic components of the deuteron one should 
require an accurate description of the bound state and, 
at the same time, the '^Si and ^Di A^A^ scattering partial 
waves. For consistency one should also require coherent 
results for the other A^A^ 5* wave, the ^Sq. The model of 
Refs. ^ fulfills the above requirements. The present 
study has been done by means of a Lippmann-Schwinger 
formulation of the RGM equations in momentum space 

In our calculation Eq. (20) of Ref. has been gener- 
alized to a coupled channel scheme. It implies a modifica- 
tion of Eq. (24) of Ref. ||2^ with additional terms which 
contain the different AA^ A^A^* couplings. The cal- 
culation of the transition potential has been simplified 
using the Born-Oppenheimer approximation. We have 
carefully checked the quality of our approximation in or- 
der to guarantee that our results will not be influenced. 
It can be seen that the on-shell properties for the low 
angular momentum partial waves obtained by means of 
the Born-Oppenheimer method are of the same quality 
(within 5% for energies below 250 MeV) as those ob- 
tained, using the same quark model Hamiltonian, but 
making use of the RGM . 

The potential yields a fairly good reproduction of the 
experimental data up to laboratory energies of 250 MeV. 
In Table | we quote the quark model parameters, that 
have been taken from Ref. [Q. As can be seen, a fine 

TABLE I: Quark-model parameters. 
mg{MeV) 313 
b(fm) 0.518 

Os 0.498 

Qch 0.0258 

m<,(fm-i) 3.538 
m,r(fm~^) 0.70 
A(fm"^) 4.30 

tuning has been done in order to consider the other non- 
nucleonic configurations. Let us finally mention that the 
quark-quark interaction used reproduces the even-parity 
baryon spectrum (A^, Af*(1440), A(1232), A(1600),..) in 
an exact Faddeev calculation |24]. 

Within the above described framework we have cal- 
culated the deuteron binding energy and wave function 
making emphasis in a simultaneous description of the 
NN scattering phase shifts. For the bound state prob- 
lem, Eq. (20) of Ref. |2^ can be discretized in momen- 
tum space and written as 

J2{H.j~ES.,)^,^0, (1) 
j 
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where we have used a simphfied notation and the indices 
i and j run not only for all the discretization points but 
also for the different channels included in the calculation. 
The nontrivial solutions are given by the zeroes of the 
Fredholm determinant 



(2) 



being the values of E that satisfy the previous equation 
the energy of the bound states. Once the energies have 
been found the wave function can be easily calculated 
solving the linear problem of Eq. (^. 

We show in Table |l| the different configurations and 
partial waves included in our calculation. In Table [II we 



TABLE II: Different channels and partial waves considereded 
in our calculation. We give on the second column the energy 
difference (MeV) with respect to the A'^A'^ system 

NN '^Si - ^Di 0.0 



iViV*(1440) 
AA ^Si 



'Si - 'Di 



501.0 



'Di - ''Di - ^Gi 586.0 



present the results obtained for the nonnucleonic prob- 
abilities and the static properties of the deuteron. In 
all cases the deuteron binding energy is correctly repro- 
duced, being Ed — —2.2246 MeV. We have shown the 
results of a calculation including only NN components, 
including NN and AA configurations and finally the 
full calculation including also A^A^*(1440) configurations. 
Among the allowed configurations, the A^A^*(1440) has 
not been usually included in the deuteron calculations 
due to the great uncertainty associated to the coupling 
constant and cutoff parameters |8[ The first result 
we would like to emphasize is the fact that the proba- 
bihty of A^A^*(1440) channels are smaller than the A A 
ones. They do not show much influence on the static 
properties of the deuteron as it seems to be case for the 
deuteron form factors However, these small compo- 
nents find support in the explanation given in the liter- 
ature to some deuteron reactions @, |[ 1). Subthreshold 
antiproton production in d — p and d — d reactions |Q , 
pd dp reactions or antiproton-deuteron annihila- 
tion at rest [Q are compatible with small percentage of 
iViV* (1440) in the deuteron wave function. 

The prediction we obtain for the 7V7V*(1440) proba- 
bihties, a larger component of the ^£>i(A^A^*(1440)) par- 
tial wave than the •^S'i(A^A^*(1440)) partial wave, agrees 
with the ordering obtained by other calculations present 
in the literature ijlll, HJ] . This can be understood if one 
takes into account that the tensor coupling, which is the 
main responsible for the presence of nonnucleonic com- 
ponents on the deuteron wave function, is much stronger 
for the ^Si{NN) -^^Di{NN*{lUO)) transition than for 
the ^Di{NN) ^^Si{NN*{lUO)) one, enhancing in this 
way the _D-wave influence with respect to the S'-wave 
component. Regarding the absolute value of the prob- 
abilities, our results are a factor ten smaller than those 
reported on Ref. plj], where an estimation of 0.17% for 



the iViV* (1440) configuration was obtained (0.06% for 
the 3 S'l and 0.11 % for the ^Di partial wave). The de- 
pendence of this result on the hypothesis made and the 
deviation from the results we obtain could be understood 
in the following way. The deuteron is calculated using the 
Reid hard-core potential. When including A^A^*(1440) 
components, the channel coupling induces an attractive 
interaction on the NN system, that needs to be sub- 
tracted out. Such a subtraction was done by reducing 
the intermediate range attraction of the central part of 
the Reid hard-core potential without modifying the ten- 
sor part, as was done in Ref. to calculate the prob- 
ability of AA components. As a consequence, in these 
type of calculations the strength of the tensor coupling 
to the A^A^*(1440) state can be enhanced by decreasing 
the intermediate range attraction in the NN channel. 
The balance between these two sources of attraction can- 
not be disentangled in a clearcut way. This is a similar 
problem to the one arising in the ^So{NN) partial wave 
when the coupling to the iVA system was included [ p6[ . 
The same attractive effect could be obtained by a cen- 
tral potential or a tensor coupling to a state with higher 
mass, being necessary other observables to discriminate 
between the two processes [26 . This seems to be the 



reason of the much bigger probability for the TV TV* (1440) 
components in Ref. jll) , that on the other hand showed 
a great dependence on the choice of the NN phcnomcno- 
logical potential. In Ref. ||2^, although the contribution 
of resonance configurations has been included to study 
the nucleon-nucleon system, there are no numerical pre- 
dictions to compare with. 

There are other estimations on the literature. The re- 
sults of Ref. [|j are only of qualitative interest. The 
pathological behavior of the transition potential to reso- 
nance states was regularized by a cutoff factor that made 
the potential too weak at small distances. In Ref. |2^ 
they study the effective numbers for different resonance 
configurations on the deuteron making use of baryon 
wave functions obtained from the diagonalization of a 
quark-quark interaction containing gluon and pion ex- 
change in a harmonic oscillator basis including up to 2huj 
excitations, and deuteron wave functions obtained from 
the Paris potential or a different quark model approach. 
They obtain an upper limit of 1% for AA components 
and 0.1% for A^A^*(1440) in agreement with the order of 
magnitude and ordering of our results. 

To summarize, we have studied the resonance struc- 
ture of the deuteron using quark-model based interac- 
tions. The A^A^*(1440) baryon resonance probabilities in 
the deuteron, being smaller than the A A ones, are in 
agreement with indirect estimations obtained from the 
analysis of processes like subthreshold antiproton pro- 
duction on d — p collisions or pd — > dp reactions. 
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TABLE III: Deuteron wave function (%) 
NN AA*(1440) AA 





















As(fm-i/2) 




95.3780 
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0.8941 
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95.1885 
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0.0148 


0.1224 
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0.1245 
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0.8941 
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